Saturation of Turbulent Drag Reduction in Dilute Polymer Solutions 
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Drag reduction by polymers in turbulent wall-bounded flows exhibits universal and non-universal 
aspects. The universal maximal mean velocity profile was explained in a recent theory. The sat- 
uration of this profile and the crossover back to the Newtonian plug are non-universal, depending 
on Reynolds number 72e, concentration of polymer c p and the degree of polymerization N p . We 
explain the mechanism of saturation stemming from the finiteness of extensibility of the polymers, 
predict its dependence on c p and N in the limit of small c p and large Re, and present the excellent 
comparison of our predictions to experiments on drag reduction by DNA. 



Introduction. The onset of turbulence in fluid flows is 
accompanied by a significant increase in the drag. This 
drag poses a real technological hindrance to the transport 
of fluids and to the navigation of ships. It is interesting 
therefore that the addition of long chained polymers to 
wall-bounded turbulent flows can result in a significant 
reduction in the drag [l| . The basic experimental knowl- 
edge of the phenomenon had been reviewed and system- 
atized by Virk 0; the increase in mean velocity profile 
as a function of the distance from the wall depends on 
the characteristics of the polymer and its concentration, 
but cannot exceed a universal asymptote known as the 
"Maximum Drag Reduction" (MDR) curve which is inde- 
pendent of the polymer's concentration or its character- 
istics. When the concentration is not large enough, the 
mean velocity profile follows the MDR for a while and 
then crosses back to a Newtonian-like profile, known as 
the "Newtonian plug" cf. Fig. ^ Recently the nature of 
the MDR and the mechanism leading to its establishment 
were explained |fj. For Newtonian flows the momentum 
flux is dominated by the so-called Reynolds stress, lead- 
ing to a logarithmic (von-Karman) dependence of the 
mean velocity on the distance from the wall With 
the addition of polymers, while momentum is produced 
at a fixed rate by the forcing, polymer stretching results 
in a suppression of the Reynolds stress (and thus of the 
momentum flux from the bulk to the wall). Accordingly 
the mean velocity in the channel must increase. It was 
shown that when the concentration of the polymers is 
large enough there exists a new logarithmic law for the 
mean velocity with a slope that fits existing numerical 
and experimental data. The law is universal, thus ex- 
plaining the MDR asymptote. It turned out that the 
polymer stretching leads to an effective viscosity which 
increases with the distance from the wall. This effec- 
tive viscosity attains a self- consistent profile, increasing 
linearly with the distance from the wall. With this pro- 
file the reduction in the momentum flux from the bulk to 
the wall overwhelms the increased dissipation that results 
from the increased viscosity. Thus the mean momentum 
increases in the bulk, and this is how drag reduction is 
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FIG. 1: Mean normalized velocity profiles as a function of 
the normalized distance from the wall during drag reduction. 
The data points from numerical simulations (green circles) |3J 
and the experimental points (open circles) p| represent the 
Newtonian results. The red data points (squares) repre- 
sent the Maximum Drag Reduction (MDR) asymptote. The 
dashed red curve represents the theory of |(| which agrees 
with the universal law Q. The arrow marks the crossover 
from the viscous linear law to the asymptotic logarithmic 
law Q. The blue filled triangles [f| and green open triangles 
represent the crossover, for intermediate concentrations 
of the polymer, from the MDR asymptote to the Newtonian 
plug. This Letter presents a theory of this crossover. 



realized. In [!j it was demonstrated by DNS that Navier- 
Stokes flows with viscosity profiles that vary linearly with 
the distance from the wall indeed show drag reduction in 
close correspondence with the phenomena seen in full vis- 
coelastic simulations. 

Theory of the saturation of drag reduction. The 

aim of this Letter is to provide a theory of the crossovers 
from the MDR to the Newtonian plug for small polymer 
concentration c p and large Tie. We show that the mech- 
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FIG. 2: Schematic mean velocity profiles. Region 1: y + < y+, 
- viscous sublayer. Region 2: y + > - logarithmic layer 
for the turbulent Newtonian flow. Region 3: y+ < y + < = 
(1 + Q) 2/n ~~ MDR asymptotic profile in the viscoelastic flow. 
Region 4: y + > y+ - Newtonian plug in the viscoelastic flow. 



anism for crossover in that limit is the saturation of the 
increase of effective viscosity: at some point there is no 
enough polymer to supply the necessary linear profile of 
effective viscosity. We will develop the theory for a chan- 
nel geometry, and then apply the results to a rotating 
disk experiment with DNA as the drag reducing agent. 
For a channel of width 2L we denote by x, y and z the 
streamwise, the wall normal and the spanwisc directions 
respectively. The only non-zero mean velocity compo- 
nent is V(y) in the streamwise direction. The traditional 
wall units are 

Ve T = L^L/is Q , y+ = yVe T /L , V+ = V/^7Z . (1) 

Here p' = —dp/dx is the pressure gradient in the stream- 
wise direction a;. 7Ze T is known as the friction Reynolds 
number, vq is the kinematic viscosity of the neat fluid. In 
terms of these variables, the predictions of the standard 
Newtonian theory and the theory of [6j are summarized 
with the help of Fig. in which the Newtonian and vis- 
coelastic velocity profiles are shown schematically in the 
reduced coordinates 0J. In the Newtonian flow there are 
two distinct regions: a viscous sublayer (region I) with 
the reduced velocity given by 

V+(y + )=y + , (2) 

and a logarithmic layer (region 2), in which velocity is 
given by Prandtl-Karman law: 

V + (y + ) = K^" 1 logy + +S N , (Prandtl-Karman law). (3) 

Here k" 1 ~ 2.29 and B N ~ 6.13 Hfl. The mean velocity 
profile in the viscoelastic case consists of three regions 



: a viscous sublayer, a logarithmic polymeric sublayer 
(region 3 in the Fig. |2J) with the slope greater then the 
Newtonian one: 

V + (y + ) = K' 1 logy+ + B v , (MDR profile), (4) 

(ft" 1 ~ 11.7, B v ~ —17 2]), and a Newtonian plug (re- 
gion 4). In the last region the velocity follows a log law 
with the Newtonian slope, but with some velocity incre- 
ment AV + : 

V+ = ft^ 1 log y + + B N + AV + . (5) 

Note that the three profiles Eqs. 10), ©, and (@J inter- 
cept at one point y + — y^ ~ n~ l ~ 11.7. 

The increment AV + which determines the amount of 
drag reduction is in turn determined by the crossover 
from the MDR to the Newtonian plug (see Fig. |2J). We 
refer to this cross over point as y"J" . To measure the qual- 
ity of drag reduction we introduce a dimensionless drag 
reduction parameter 

y+ 

Q=^L-1. (6) 

The velocity increment AV + is related to this parameter 
as follows 

AV+ = (ft^ 1 - ft; 1 ) log (y+/y+) = alog(l + Q) . (7) 
Here a = k~ x — k" 1 ~ 9.4. The Newtonian flow is then 

V N 

a limiting case of the viscoelastic flow corresponding to 
Q = 0. 

The crossover point y^ is non-universal, depending on 
Tie, the number of polymers per unit volume c p , the 
chemical nature of the polymer, etc. According to the 
theory of drag reduction [6j, the total viscosity of the 
fluid ^tot(y + ) = vq + v p {y + ) [where v p {y + ) is the poly- 
meric contribution to the viscosity] is linear in y + in the 
MDR region: 

vtot(y + ) = ^oy + /?y+ , y+ <y+ <y+ . (8) 

When the concentration of polymers is small and Tie is 
large enough, the crossover to the Newtonian plug at y+ 
occurs when the polymer stretching can no longer provide 
the necessary increase of the total fluid viscosity. In other 
words, in that limit the crossover is due to the finite 
extensibility of the polymer molecules. Obviously, the 
polymeric viscosity can not be greater than v v max which 
is the viscosity of the fully stretched polymers. Thus 
the total viscosity is limited by v$ + i/ pmM . Equating 
and vtotiVy ) gives us the crossover position 

vt = ("o + "p™)/"" ■ ( 9 ) 

It follows from Eq. © that the drag reduction parameter 
is determined very simply by 

Q = ^pmra/fo , c p small, Tie large . (10) 
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At this point we need to relate the maximum poly- 
meric viscosity v pmax to the polymer properties. To this 
aim we estimate the energy dissipation due to a single, 
fully stretched, polymer molecule. In a reference frame 
co-moving with the polymer's center of mass the fluid ve- 
locity can be estimated as u ~ rVu (the polymer's center 
of the mass moves with the fluid velocity due to negli- 
gible inertia of the molecule) . The friction force exerted 
on the z-th monomer is estimated using Stokes law, 



(11) 



where a is an effective hydrodynamic radius of one 
monomer (depending on the chemical composition), and 
Ti is the distance of the i-th monomer from the center 
of the mass. In a fully stretched state ~ ai (the 
monomers are aligned along a line). The energy dissi- 
pation rate (per unit volume) is equal to the work per- 
formed by the external flow 
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^pmax = v a 3 c p N^ . (13) 



Finally, the drag reduction parameter Q is given by 



Q 



a 3 c 



p jv p 



c p small, Tie large 



(14) 



This is the central theoretical results of this Letter, re- 
lating the concentration c p and degree of polymerization 
iV p to the increment in mean velocity AV + via Eq. J71. 

Testing and explaining a DNA experiment 11]. 
A particularly interesting experiment suitable for test- 
ing our prediction was described in [TT|| . Here turbu- 
lence was produced in a rotating disk apparatus, with A- 
DNA molecules used to reduce the drag. The Reynolds 
number was relatively high (the results below pertain to 
Tie « 1.2 x 10 6 ) and the initial concentrations of DNA 
relatively low (results employed below pertain to 2.70 
and 1.35 wppm). During the experiment DNA degrades; 
fortunately the degradation is very predictable: double 
stranded molecules with 48 502 bp in size degrade to dou- 
ble stranded molecules with 23 100 bp. Thus invariably 
the length N p reduces by a factor of approximately 2, 
and the concentration c p increases by a factor of 2. The 
experiment followed the drag reduction efficacy measured 
in terms of the percentage drag reduction defined by 



FIG. 3: %DR in a rotating disk experiment with A-DNA as the 
drag reducing polymer. Note that the %DR is proportional 
to c p . When the length iVp reduces by a factor of 2 and, 
(\2) simultaneously, c p increases by factor of 2, the %DR reduces 
by a factor of 4. 

without and with polymers, respectively. The main ex- 
perimental results which are of interest to us are summa- 
rized in Fig. We see from the experiment that both 
initially (with un-dergraded DNA) and finally (with de- 
graded DNA) the %DR is proportional to c p . Upon de- 
grading, which amounts to decreasing the length N p by a 
factor of approximately 2 and, simultaneously increasing 
c p by factor of 2, %DR decreases by a factor 4. 

The flow geometry is rather complicated: with a rotat- 
ing disk the linear velocity depends on the radius, and 
the local Tie is a function of the radius. The drag re- 
duction occurs however in a relatively small near-wall 
region, where the flow can be considered as a flow near 
the flat plate. Thus, we consider an equivalent channel 
flow - with the same Tie and a half width L of the order 
of height/radius of the cylinder. In this plane geometry 
the torques in (|15fl should be replaced by the pressure 
gradients p' N v : 



T -T 
n 2 v 
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x 100. 



(15) 



where T N and T v are the torques needed to maintain 
the disk to rotate at a particular Reynolds number Tie 



= ^ ■ Fv x 100 



(16) 



In order to relate %DR with the drag reduction parameter 
Q, we re- write Eq. (|SJ in natural units 



K N J log [v\/Vl/vg 
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With constant Tie the centerline velocity Vq = V(L) is 
kept fixed: 

%s M =7& r /e -i] og ft BT + B +AV+] . (18) 
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FIG. 4: Drag reduction efficiency %DR as function of the 
drag reducing parameter Q for different Reynolds numbers 
Tie: 1.2 x 10 5 , 1.2 x 10 6 , and 1.2 x 10 7 (from top to bottom). 



This equation implicitly determines the pressure gradient 
and therefore the %DR as a function of Q and Tie. The 
set of Eqs. (7J and (|18fl is readily solved numerically, and 
the solution for three different values of lie is shown in 
Fig. El The middle curve corresponds to Tie — 1.2 x 10 6 , 
which coincides with the experimental conditions 
One sees, however, that the dependence of %DR on the 
Tie is rather weak. 

One important consequence of the solutions shown in 
Fig. His that for small Q (actually for Q < 0.5 or %DR < 
20), 93DR is approximately a linear function of Q. The 



experiments 1 1 1| lie entirely within this linear regime, in 
which we can linearize Eq. (|18|) and find an approximate 
solution for the %DR: 



%DR: 



2aQ 



- 1 log(e7&°) + 



x 100 



(19) 



Here 7<fe° is the friction Reynolds number for the Newto- 
nian flow, i.e. the solution of Eq. (|18fl for AV + = 0. 

It is interesting to note, that while the %DR depends on 
the Reynolds number, the ratio of different %DR's does 
not [to 0(Q)\. 
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This result, together with Eq. (jl3|) . rationalizes com- 
pletely the experimental finding of 11] summarized in 
Fig. |21 During the DNA degradation, the concentration 
of polymers increases by a factor of 2, while the num- 
ber of monomers iVp decreases by the same factor. This 



means that %DR should decrease by a factor of 4, as is 
indeed the case. 



Summary. We have presented a theory of the 
crossover from the universal MDR mean velocity profiles 
to the Newtonian plug. We have connected the theory 
to measured percentages of drag reduction, and tested 
the theory against experiments in which A-DNA is used 
as the drag-reducing agent. The experimental results 
pertain to high Tie and small c p , where we can assert 
that the crossover results from exhausting the stretching 
of the polymers such that the maximal available viscosity 
is achieved. In the linear regime that pertains to this 
experiment the degradation has a maximal effect on the 
quality of drag reduction Q, leading to the precise factor 
of 4 in the results shown in Fig. |3] Larger values of the 
concentration of DNA will exceed the linear regime as is 
predicted by Fig. 01 then the degradation is expected to 
have a smaller influence on the drag reduction efficacy. 
It is worthwhile to test the predictions of this theory also 
in the nonlinear regime. 
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